Exact Unitary Transformation of the 
One-Dimensional Periodic Anderson Model 
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An effective hamiltonian is derived exactly for the one-dimensional periodic Anderson model via 
a canonical transformation. The canonical transformation has been calculated up to infinite order, 
thus an exact transformation was performed in the strict mathematical sense. We also discuss briefiy 
the impact of the obtained result on understanding the magnetic properties of several Kondo lattice 
compounds. 
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The periodic Anderson model (PAM) is believed to 
contain the essential physics to describe the low tem- 
perature magnetic, superconducting or semiconducting 
properties of the heavy fermion materials and yet the 
understanding of it is still very limited. The model is 
exactly solvable via Bethe Ansatz only in the single im- 
purity case m and even in one-dimension (ID) it's Kondo 
lattice limit only allows an exact solution via bosoniza- 
tion 0. 

To begin understanding the role that PAM plays in 
determining the properties of heavy fermion systems we 
must gain a better understanding of the effective inter- 
actions and local correlation present in the model. We 
achieved this by applying a canonical transformation to 
PAM, similar to that used by Schrieffer and Wolff 0, 
which has been originally constructed in order to estab- 
lish the connection between the single impurity Anderson 
model and the single impurity Kondo model. 

The Schrieffer- Wolff transformation Q was performed 
only up to first order, which restricts the validity of the 
transformation to a very limited range of the original pa- 
rameters. In this Letter we focus on the ID case, where 
we developed a method that allows us to derive recursive 
equations for arbitrary order of the canonical transfor- 
mation, and enabled us to sum up the transformation 
to infinite order. Thus the canonical transformation is 
exact, and it is not suppose to eliminate any irrelevant 
degrees of freedom. After the transformation is carried 
out a new Hamiltonian is derived and the exact micro- 
scopic expressions for the effective interaction constants 
obtained. 

From these coupling constants, we will analyse in more 
detail the value of the effective magnetic coupling be- 
tween the impurities (the so-called Kondo coupling) to 
gain a deeper insight into the magnetic property of the 
ID PAM. This problem has generated considerable inter- 
est since Moller and Wolfe [5| have shown, using a slave 
boson approximation, that both ferro- and antiferro- 
magnetism may exists in the ID PAM. Recent numerical 
calculations [6| have also shown the existence of ferro- 
magnetic phase in the PAM for a wide range of dop- 



ing ratios between quarter-filling and half-filling. On 
the contrary, the ground state of the symmetric ID 
PAM is antiferro-magnetic in the half-filled case and 
favours antiferro-magnetic interaction in the vicinity of 
quarter filling 0. In addition to this, there are many 
cerium mixed valence materials 01, which exhibit ferro- 
magnetism in the presence of strong Kondo- like behavior. 
The exact Kondo coupling, as calculated hereafter, will 
shed light on this complex magnetic behavior of ID PAM. 

1. Canonical Transformation. We write the standard 
PAM Hamiltonian in IB as H = Hq + Hy, with 

Ho = -t^(cj+i + h.c.) + AiXl'^Uc^,^ 



and Hv = VJ2iAfLc.,a + h-c). Here c],^ (c,^^) create 
(annihilate) conduction electrons with spin a at lattice 
site i, fl^ ifia) create (annihilate) /-orbital impurity 

electrons (n{^ ~ ficrfia)^ ^^'^ i^f) energies 
of the conduction (/-orbital) electrons, respectively, t 
is the nearest- neighbour hopping matrix element, U the 
on-site Coulomb repulsion of the /-electrons, and V is 
the strength of the on-site hybridization matrix element 
between electrons in the /-orbitals and the conduction 
band. 

Its unitary transform is simply written as: 



H 



Ho + [S,Hv]/2 - 



[S, [S,Hv]]/3 



(2) 



where S is determined by the condition Hy + [S, Ho] = 0. 
One can easily verify that 



S = 



ViA + Zfl_J^_^ if^c,^, - c\ j,J , (3) 



with A = l/{2t-n + ef ) and Z = l/{2t~fi + ef + U)-A 
satisfies this condition. In Hq of Eq. Q a continuum rep- 
resentation (similar to the field theoretical bosonization 
scheme Q) for the kinetic energy has been used. 
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TABLE L The coefficients J, P, J, K and M (a"'^ = 71°=^") 
of the transformed Hamihonian in the first, third and fifth 
order, corresponding to the first, second and third row. The 
G coefficient has only one term in each order: 2 A in the first, 
-8^^ in the second and 32^4^ in the fifth order, respectively. 



J 


P 


/ 


if 


M 


2a°i 




2a«i 








-32a2^ 


-8a2i 


-28a^^ 


-4a2i 




-32ai^ 




-36a^^ 


+4ai2 


+8aO^ 




+2a°3 








512a*^ 


224a'*i 


336a*^ 


1760*1 


-352a32 


+ 1024a-^2 


+448a^2 


+848a^^ 




-528a2^ 


+ 10240^^ 


+328a2^ 


+936a2^ 


+88a23 


-368a" 


+5120^** 


+ 104a" 


+520a" 


-8a" 


-96a°5 


+128a''5 


-4a«5 


+128a'''^ 







Having found S*, it is straightforward to determine the 
first [n — 1), third (n — 3) and fifth [n = 5) order 
terms of the transformed Hamiltonian which depend on 
the following commutation, 



n times 



[[S,Hv]]n = [SAS,[S.---AS,Hv]. 



(4) 



This commutation was found to have the following form 
over the first three odd values of n, 



E 



'■1,17 



t, — a t,a 



(5) 



the first three odd n are summarized in Table Besides 
the terms which renormalize the starting Hamiltonian, 
we obtained three new effective interactions: J, the well- 
known Kondo coupling; P, a Josephson type two particle 
intersite tunnelling; K, an effective on-site Coulomb re- 
pulsion for the conduction electrons; and a higher order 
term, M . From these, K and M are absent in the first or- 
der Schrieffer- Wolff result, as they depict a higher order 
coupling between the impurity and conduction electrons. 
This coupling becomes stronger as U increases, but only 
to a certain extent before the coupling decreases. 

It is not by incident that the first three odd order com- 
mutation Eq. is given by a close form Eq. (0). In fact, 
we will prove in the following using induction that the 
form of the n odd order commutation Eq. Q is given by 
Eq. jSJ in general. 

Other than the form of the commutation, one can also 
find a pattern in the coefficients of the common terms 
among the first, third and fifth orders, from TableD This 
pattern has been verified true up to eleventh order and 
can be proven valid for any order by the same induction. 

2. Proof by Induction. The n = 1 case is the well 
known Schrieffer and Wolff result j3| which represents 
the first rows of Table By visual inspection, it can be 
observed that there is a pattern in the coefficients of the 
common terms among the first, third and fifth orders. 
This pattern can be proven to exist for any order by 
induction, if the same pattern remains after commuting 
Eq. lO with 5 twice. 

To do this, two different indices are introduced to dif- 
ferentiate the order of the commutation n from the re- 



currence of the coefficients Jm, P„ 



over odd orders. 



The values of the coefficients J, P, G and /, M for 



The mapping of the two sequences can be written as 
n = 2m -|- 1 for odd order n. Assuming that Eq. © 
is true for any m, we calculate its commutation with S 
for [[S*, iJy]]„+i, which yields: 



E 



-2VJm{A + Zni_^) (nl_^ - <_J - 2VP^m{A + ^<_J - <_ J - 2VGra{A + Zn^) 



-2VM^{A + Z)nlnl_^ - 2VIm{A + + 2VKUA + 



r 



f + c 

i,aJi,{T ' J 1,(7 I, a 



(6) 



-Gm4:V^{2AZ + Z^) - Irn4:V'^{A + Zf 

-K„AV^A^ , 

Pm+l = -Xm^V'^iA^ + AZ) - Pra2V''{{A + zf + A^) 

-{I.^ + K,n)2V\A + Z)A, 
J„+i = ^J,^W\iA + Zf + A^)-P,n8V^{A^+AZ) W = -J„,W^{A + Zf - P^W^{A^ + AZ) 



In the next step of the proof, calculating the commu- 
tator of Eq. ||HJ) with S again reveals that the obtained 
[[S, Hv]]n+2 has exactly the same form as [[S, Hv]]m with 
coefficients: 
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-G^AV^{2AZ + Z^) - I,nAV^{A + Zf , 
Km+i = ~J,n^V^A^ - P„AV^{A^ + AZ) - K„,AV^A^ , 
Mm+i = +K,r,AV^{2AZ + Z^) - M„AV^{A + Zf 

+Jr,AV^{'2AZ + Z^) , 
G,„+i = -G,„4y2^^ (7) 

This proves that Eq. (j^l is true for any m. Thus, using 
mathematical induction, we conclude that the form of the 
nth commutation of Hy with S is closed and is always 
given by Eq. 

4- Evaluating the Coefficients. The recursive equations 
(0 can be solved simultaneously to give the odd order 
coefficients of the transformed Hamiltonian, Jm, Pm, Inn 
Km and Mm- 

These recursive equations can be summarized into an 
matrix form in which 



^ Jm+1 


] 






\ 






-/3^\ 


^Pm+l 






P 











Im+1 
















\ Km+l 


) 




\Km 


) 




\ 


1 



(8) 



where a = A + Z and f] = A. M is a matrix given by: 



a 



M = 



2a(3 a^ + P^ a/3 a/3 



a/3 
a(3 



a 






/32 / 



(9) 



These equations form a set of simultaneous recursive 
equations with the first order, ie, the Schrieffer and Wolff 
result as the initial condition, and their solution is: 



Jrr 
Prr 

G„ 
K„ 



= {-2V^Y\a^ + (3^r~\2'^"'+^a[3 

= (-2T/2)m^2/32)"2/3y2 ^ 

= {(-2T/2)'"(a^+/32)"-i[4™/3(a-/3) 

-a(a + 13)] + (-2F^)"(2/32)™}2/3y2 ^ 
= (-23v^2)"(a2+/32)"-i2a2(a-/3)F^ 
-2F')"(a2 + /32)™-i2a/3(a + p)V^ 



-21/2)'"(2/32)™2/3y2 ^ 

-23T/2)™(a2 + /32)™-i2(a - p){a' - p^)V^ 
-2V^na'^ + p^)'"~Ha(3{a + p)V'^ 
+ (-4aV2)™2al/2 + (-4/3V2)"2/3y2 . 

The even order coefficients can then be deduced from 
the odd order coefficients, by applying another commu- 
tation. These coefficients also have a form, as one would 
expect, but the pattern is quite different to that of odd 
order. The (n + 1)*'* commutation result can be summa- 
rized in the form: 

[[S, Hv]]n+1 ^(^m + Smn{_„ + Tmn1_„ 
i,(T 



(10) 



where 



Rm = -i-W^P^)"'4P^V^ , 
Sm = -[-2^V^{a^ + /3^)r8a{a- P)V^ 
-[-21/2(^2 _^ ^2)p2/3(a + /3)V^ 

Tm = [-2^V^{a^ + (3^)r8(3ia- p)V^ 
-[-21/2(^2 ^ ^2)p2a(a + P)V^ 

Qm = [-2'v\a^ + f3^)r8ia- Pfv"" 
+ [-21/2(^2 + /32)]"2(a + /3)2y3 

-(-4T/2a2)'"4aV^ - {-AV^ P^)"Ap^V^ . 

Accordingly, we have obtained the result of the commu- 
tation of Hy with S to any order, and hence have all the 
information needed to re-build the Hamiltonian after the 
transformation. 

4- The Transformed Hamiltonian. Using the general 
expression for the nth and the (n -I- 1)*'' commutation 
of S with Hy, the infinite order transformation can be 
calculated. The transformed Hamiltonian comprises Hq, 
the sum of the odd order commutations of S with Hy 
and the sum of the even order commutations: H — Hq + 
Hodd+Heven, wi^ere Hodd = I]m=o[l/(2"^+l)'^ 1/(2?^+ 
2y.][[S,Hy]]2m+i is: 

Hodd = ^J{c\aCi,-<Tfl-afi,a-nl„'nl_a) 

~^-^i^i,(T^i,-(Tfi,(Tfi,-(T ~^ fi,crfi,-(T^i,cr^i,-(T) 

(11) 



/ / 



J, P, G, /, K and M are the summation of the cor- 
responding Jm, Pm, Gm, I-m, Km and A/,„ ovcr infi- 
nite number of to. If we define a = A + Z, /3 = A, 
e = ./WHc^'TW), 9p = 2VP, = 2ya and F(x) = 
sinx/a; -I- (cos a; — l)/a;2 then the exact values of the cou- 
pling constants from Eq. are: 



J 2{a~ l3)V'^F{2e) , 

P = 2af3{a ^ P)V^F{2e) - {a - (3)V^ 

G = 2f3V^F{e0) , 



a2 -F/32 



m 



I = 2a 



K = 2p'V'^-^F{2e)-2a(3V'^±^F{e)+G 
a'^ + p'^ a^ -f /3^ 

a ~ (3 „.„„^ „ ^,,2 OL + 13 



M = -2(a-^) 



a2+/32 
a2-/32 



F{20) + 2apV^^^-^F{e) - G , 
a'^ + /3^ 

V^F(2e) ~ Aal3^^^V^F{e) 



a2 /32 ^ ^ ' a2 + ^2 
+2a^^2^(6'„) + 2(3V^F{ep) . 
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Similarly, the even order Hamiltonian Heven = 
E"=o[l/(2m + 2)! - 1/(2to + 3y.][[S, Hv]hm+2 can be 
evaluated by substitution, 

i.<7 

x(cU4. + /LO, (12) 

where i?, S*, T and Q are the summation of the corre- 
sponding Rrrn Srn, Qm and Tm over infinite number of 
m. Using the same notations for a, /3, 9, Op, 6a as above, 
and F'{x) ~ sinx/x^ — cosx/x^, we obtain: 

R = -Al3'^V^F'{9p) , 

S = -8{a- /3)aV^F'{20) ^2{a + (3)(3V^F'{e) - R, 
T = %l3{a- l3)V^F'{2e) -2a{a + l3)V^F'{e) - R , 
Q = 8(a- (3fV^F'{29) + 2{a + pfV^F'{9) 
-4:(3V^F\9fj) - ApV^F'{9c) • 

In the symmetric case (e/ = — C//2), one gets Z — 
-2A, 9 = 9p ^ -9a ^ 2AV, a = -A and /3 = A, both 
the odd and even order Hamiltonian coefficients simplify 
considerably. 

5. Discussions. We have shown that an unitary trans- 
formation of the ID PAM can be calculated to infinite 
order. The obtained series for every coefficient gener- 
ated by the transformation are convergent. Even though 
the transformed Hamiltonian appears to be more compli- 
cated than the original one, it gives the exact expression 
of the effective interactions which are valid for any U and 
V. These represent a vital source of infromation due to 
the lack of similar, exactly soluble theories of PAM. 

The Kondo coupling J from our exact result has a 
behaviour different to that of the well known Schrieffer- 
Wolff result. This difference is attributed to the higher 
order terms of the coupling coefficient J in the trans- 
formation. As depicted in Fig. ^ these higher order 
terms seem to have an effect of suppressing the anti- 
ferromagnetic coupling between the localized electrons 
for large values oiV/U. As V/U increases, J can change 
sign, indicating a possible change of the magnetic cou- 
pling state which has been observed in numerical ap- 
proaches 0, . The value of J is also crucial to deter- 
mine ^3 the Kondo (oc exp — 1/| Jp|) and Neel (oc | Jpp) 
temperatures, where p is the density of states for the con- 
duction band at the Fermi level. In the Kondo regime, 
see Fig. ^ J passes through a maximum, as a function 
of U and V. Accordingly, the Kondo and Neel temper- 
atures will show a similar behavior. Such behavior has 
been observed experimentally with increasing pressure 
(i.e., V) in several cerium Kondo compounds |9|. Hence, 
this exact result will also shed new light in understand- 
ing the competition between the Kondo effect and the 



Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in 
ID and revisit the Doniach diagram 

HI 

0.9 - - 

0.8 - y-' - 

0.7 - ^ ^ 

0.6 - y^ 



5 0.5 - 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

8V/U 



FIG. 1: The spin coupling constant J up to different order 
of transformation. HI is the Schrieffer Wolff result, H1-I-H3 
is the result up to the third order of the transformation, and 
similarly for H1-I-H3-I-H5. Solid line shows the result up to 
infinite order. 



Finally, one may notice that the hybridization term 
is still present in the final result of the transformation, 
see Eq. \\2\ . However, the strength of the hybridiza- 
tion is substantially reduced to the third order of V as 
--V{{w^2VA)l2VA- co%2VA) x -^A^V^ for small 
V. It is in fact possible to remove this hybridization 
term completely, but the transformation involved will be 
slightly different [ij. 
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